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Resumen 


Se describe la calidad química de los recursos hídricos tanto subterráneos como 
superficiales en la cuenca del río Duero. El agua subterránea se analizó mediante 
un estudio hidrogeoquímico en 97 aprovechamientos (20 manantiales y 76 pozos) 
para la época de estiaje. Se aplicó el Índice de Calidad del Agua (ICA; NSFWQI) 
en el agua superficial para 35 sitios, a partir de ocho parámetros, dividiéndolos 
en manantiales y el cauce. Se compararon diferentes parámetros con información 
previa mediante un análisis de varianza. La calidad del agua subterránea en términos 
generales resultó buena y está asociada con el tipo de rocas y geología de la cuenca. 
Predomina la clase de agua C1-S1, baja en salinidad y sodio, pudiéndose utilizar 
en la mayor parte de los cultivos y cualquier tipo de suelos sin desarrollar peligro 
por salinidad y sodicidad. Los valores ICA obtenidos, uso en agua potable, ubican 
a la mayoría de los manantiales como contaminados, con excepción del manantial 
de Carapan (condición excelente con 92 puntos). En el cauce, uso en agricultura, el 
75% de los sitios resulta como contaminado y el resto en levemente contaminado. 
La comparación de los nitratos en 20 años de diferencia fue estadísticamente 
significativo (F = 15.73, p = 0.0001), lo que representa un aumento importante de 
los nutrientes en el río. Por último, el creciente deterioro de la calidad del agua del 
río, a pesar de un balance hidrológico positivo, promueve la extracción del agua 
subterránea, llevando a un proceso de sobreexplotación del manto acuífero. 


Palabras clave: análisis hidrogeoquímico, índice de calidad del agua, contaminación, 
sobreexplotación. 


Introducción 


En México, el Plan Nacional de Desarrollo 2006- 
2012 incluye el Programa Nacional Hídrico, 
en el cual se señalan los pasos para un manejo 
adecuado del recurso de manera comunitaria 
(Presidencia de la República, 2007). Sus ob- 


jetivos principales son prevenir y reducir la 
contaminación del agua; fomentar su uso 
adecuado; proteger ríos, lagos y lagunas; y 
reducir los efectos de inundaciones y sequías 
(Conagua, 2008). 

Para lograrlo, la Comisión Nacional del Agua 
(Conagua) ha promovido la realización de lo 
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que denomina proyectos emblemáticos para el 
saneamiento integral de cuencas, considerados 
de alto beneficio social y que pueden ser 
concluidos en cinco años. Con ello se busca 
dar respuesta más efectiva a la solución de la 
problemática socioambiental de cuenca, con 
base en un consenso de usuarios. 

Esto responde a que en la actualidad podría 
establecerse que prácticamente la totalidad 
de las cuencas hidrológicas de México se 
encuentran en mayor o menor medida con 
verdaderos problemas de contaminación. 
Ejemplo de ello es la cuenca del río Lerma, 
considerada la más contaminada de México 
(Pacheco-Vega, 2007). 

La cuenca del río Duero, localizada dentro 
de la zona hidrológica río Lerma-Chapala no 
escapa a esta situación. Problemas relacionados 
con la pérdida de suelo, disminución de la 
cubierta vegetal, tala clandestina de bosques, 
contaminación de suelos y acuíferos por 
basura, descarga de aguas negras directas al 
río, y desaparición de fauna silvestre, son el 
común denominador (Conagua, 2009). 

Dado el deterioro continuo en la calidad 
del agua superficial de la cuenca del río Duero, 
el aprovechamiento del agua subterránea 
en la misma se ha incrementado de forma 
considerable. El Registro Público de Derechos 
de Agua (Repda) de la Conagua registra la 
existencia de un total de 150 pozos para el 
año 2000, pasando en el 2012 a 805 pozos, 
con volúmenes totales de extracción anual 
de aproximadamente 120 Mi. En la cuenca, 
los acuíferos presentan condiciones de libre 
a semiconfinado, lo cual la hace altamente 
vulnerable a actividades antropogénicas. 

En atención a dicha problemática, es que se 
desarrolla el presente estudio, con el objetivo 
principal de evaluar la actual problemática del 
agua subterránea y superficial en la cuenca, en 
términos de su calidad química, que impacta 
una importante actividad económica en el 
ámbito estatal en el sector agrícola, con el 
cultivo de frutillas, hortalizas y granos, que 
incluyen productos de exportación, como la 
fresa. 


Metodología 


Área de estudio. La cuenca del río Duero se 
ubica al noroeste del estado de Michoacán, 
en las coordenadas 19” 40” y 20* 15” latitud 
norte, y 101% 45 y 102” 45 longitud oeste. 
Pertenece a la Jurisdicción Administrativa 
núm. VIII Lerma-Santiago-Pacífico, dentro 
de la Región Hidrológica núm. 12 Lerma- 
Santiago, específicamente en la denominada 
Subregión Hidrológica Bajo Lerma. Tiene una 
superficie de aportación de 2 531.3 km? con un 
escurrimiento promedio anual de 250 millones 
de metros cúbicos. Tiene forma alongada y 
corre en dirección sureste-noroeste, cruzando 
cuatro microcuencas: la cañada de los Once 
Pueblos, los valles de Guadalupe y Zamora, 
para, finalmente llegar hasta la Ciénaga de 
Chapala (figura 1). 

Análisis de la calidad del agua subterránea. Con 
el fin de determinar el esquema de distribución 
de la calidad química del agua subterránea 
para uso agrícola en la cuenca del río Duero, 
se realizó un muestreo hidrogeoquímico en 
96 distintos tipos de aprovechamientos (20 
manantiales y 76 pozos) durante el periodo 
abril-mayo de 2009. 

Durante el muestreo se obtuvieron datos en 
campo de temperatura, conductividad eléctrica, 
sólidos totales disueltos, y pH, utilizando un 
conductivímetro portátil de medición, marca 
HANNA Instruments, modelo HI 9635. 

En la determinación de los principales ele- 
mentos químicos (Ca*, Mg*, Na”, K* HCO",, 
Cr y SO?”,) necesarios para la obtención 
de las distintas facies y clases de agua, se 
siguieron patrones establecidos en el manual 
de Standard Methods for Examination of Water € 
Wastewater (APHA, 1998) y fueron realizados 
en el Laboratorio de Suelos y Aguas del 
Centro Interdisciplinario de Investigación 
para el Desarrollo Integral Regional, Unidad 
Michoacán (CHDIR-IPN-Michoacán). 

Con el Diagrama de Piper se definieron las 
distintas facies químicas de agua prevalecientes 
en la cuenca. Las aguas geoquímicamente 
similares quedan agrupadas en áreas bien 
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Figura 1. Mapa de localización de la cuenca del río Duero. 


definidas. Las concentraciones de los iones 
dadas en meq/L, se expresan como porcentaje 
con respecto a la suma de aniones y cationes, 
respectivamente. 

Para definir el peligro de salinización o 
sodificación del suelo a partir de parámetros 
medidos en el agua de riego, se utilizó el 
esquema propuesto por el U.S. Soild Salinity 
Laboratory Staff, en 1954, conocido como 
Normas de Riverside. Este sistema se basa en 
la medida de la conductividad eléctrica del 
agua para determinar el riesgo de salinización 
del suelo y en el cálculo de la Relación de 
Adsorción de Sodio (RAS), para definir el riesgo 
de sodificación o alcalinización, mediante la 
fórmula: 


Na'* 
Ca'*+Mg'* 
2 


RAS= 


Donde las concentraciones de Na', Ca** y 
Mg” están dadas en meq/L. 

Se definen cuatro clases de riesgo de 
salinización y otras cuatro de riesgo de 
sodificación, resultando en total 16 clases 
de agua (C1-S1, C1-S2, etc.) cada una de las 
cuales tiene diferentes condicionantes para ser 
utilizada en regadío. 

Se elaboraron los mapas de distribución 
de isovalores ¡para los parámetros de 
conductividad eléctrica (CE), sólidos totales 


disueltos (STD), temperatura y relación de 
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adsorción de sodio (RAS), aplicando en todos 
los casos el método de interpolación Kriging. 

Análisis de la calidad del agua superficial. Se 
muestrearon 35 sitios en el río Duero, incluyen- 
do el cauce principal, afluentes secundarios y 
manantiales. Dentro de las variables ambien- 
tales se midieron la temperatura (“C), pH, 
oxígeno disuelto (ppm), % de saturación 
de oxígeno y nitratos (ppm) por medio de 
un multisensor (Hydrolab DataSonde 4X). 
Además, se recolectaron 2.5 L de muestra en el 
primer metro de profundidad para determinar 
la alcalinidad, demanda bioquímica de oxígeno 
(DBO,; en ppm), demanda química de oxígeno 
(DQO; en ppm), sólidos totales suspendidos 
(ppm) y fosfatos (APHA, 1998). Por último, 
se estimaron las bacterias coliformes fecales, 
mediante lo establecido en la NOM-112- 
S5A1-1994, y según Vanderzant y Spittstoesser 
(1992), empleando la técnica del número más 
probable (NMP). 

Para la determinación del índice de calidad 
del agua (ICA), se incluyó la medición de ocho 
variables (temperatura, pH, oxígeno disuelto, 
coliformes fecales, demanda bioquímica de 
oxígeno, fosfatos, nitratos y sólidos totales 
suspendidos). Se siguió el método del Water 
Quality Index desarrollado por la National 
Sanitation Foundation (WOINSF), en donde 
la evaluación se realiza a partir de técnicas 
multiplicativas y ponderadas. Esto asigna un 
peso específico a cada una de las variables, 
según la importancia que tenga para la vida 
acuática y el ambiente. Se aplicó la ecuación 
propuesta por Brown et al. (1973): 


n 
ICA=| [y;" 
i=1 


Donde n son las ocho variables; w, el peso 
específico que se asigna a cada parámetro (1), 
que se pondera entre 0 y 1, de tal forma que 
se cumpla que la sumatoria sea igual a uno; q, 
la calidad del parámetro (i), en función de su 
concentración, y cuya calificación oscila entre O 
y 100; PI representa la operación multiplicativa 
de las variables Q elevadas a la W. 


Las gráficas de sensibilidad, en donde, como 
función de la concentración de la variable, se 
lee la calificación de la calidad y se obtienen 
de las referencias originales. En el caso de la 
interpretación se dividieron los sitios entre 
manantiales, donde se aplicó el criterio para 
uso de agua potable, y el cauce del río con el 
criterio de uso agrícola. 

El valor ICA obtenido se acomodó en una 
escala entre 0 y 100, que califica la calidad, a 
partir del cual y en función del uso del agua, se 
puede estimar el nivel de contaminación. Esto 
permite interpretar y comunicar fácilmente los 
resultados en una escala sencilla, agrupando 
la calidad del agua de los ecosistemas acuá- 
ticos en seis intervalos: (E) Excelente; (A) 
Aceptable; (LC) Levemente Contaminada; (C) 
Contaminada; (FC) Fuertemente Contaminada, 
y (EC) Excesivamente Contaminada. 

Dado que se cuenta con información previa 
de algunas de las variables (López et al., 1990), 
se hizo una comparación de la información 
de manera gráfica y se realizó un análisis de 
varianza de una vía (ANOVA) para determinar 
si hay diferencias significativas estadísticas. 
Se determinaron las diferencias mediante la 
prueba de Tukey, con un intervalo de confianza 
del 95% (P < 0.05), por medio del programa 
SPSS Statistics, ver. 17.0 (2008). 


Resultados 
Agua subterránea 


La figura 2 muestra la localización dentro de 
la cuenca de los 76 pozos y 20 manantiales 
muestreados. 

Los resultados de los análisis químicos 
realizados en manantiales y pozos se muestran 
en el cuadro 1 y 2, respectivamente. 

Para tener una idea preliminar de la 
calidad química del agua para uso agrícola 
se obtuvieron mapas de isovalores para los 
parámetros de conductividad eléctrica (CE), 
sólidos totales disueltos (STD) y temperatura 
(T). La figura 3 muestra la configuración de los 
mismos en la cuenca. 
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Figura 2. Localización dentro de la cuenca de los distintos aprovechamientos de aguas subterráneas muestreados: 
a) pozos, b) manantiales. 


Cuadro 1. Resultados de análisis físico-químicos a muestras de agua de manantiales. 


A E O OS K Mg Ca Cr | HCO,; | CO; | sos 
(us) (ppm) (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) 

Carapan 131 93 18.2 23.3 17.9 72 8.2 9.2 110.5 0 2.1 
Tanaquillo 148.6 105.8 AS DNS 19.3 Yen 10.2 11.3 des 0 2.6 
Chilchota 191 133 19.4 25.4 23.4 8 9.7 11.3 114.9 0 5.9 
Nogales 1932 96.5 18.1 pol US 7.8 9.6 8.5 116.1 0 222 
Etúcuaro 153 107.5 19.4 25.1 21.3 LO 9.9 13.5 112.8 0 2.8 
Telonzo 178 105.8 237% 24.9 27.4 Ol) 10.5 1519 159211 0 24) 
Tangaman 181 126.7 22.9 26.8 25.6 9 10.7 12,3 128.9 0 6.8 
Estancia 190 199:2 22218) 202 25.6 8.2 10.2 10.1 128.3 0 2.4 
Orandino 139 98 21.9 23.5 29.8 8.2 9.8 11.5 125.7 0 2.3 
Verduzco 155 108 2 DAD) 28.8 8.2 10.2 12.8 128.5 0 25) 
Platanal 190 133.2 22.5 28.5 27.9 7.3 10.8 12.8 129.1 0 2.3 
Camécuaro 243 170 22.4 30.8 28.7 8.1 10.6 20.9 126.1 0 3.4 
Cupatziro 101 78 18.3 28.8 26.7 6.4 11,1 12.2 129.6 0 2.2 
San José 85 60 18.2 ¿LS 18.4 6.4 8.9 7.8 103.3 0 4.9 
Antonio 253 179 22.5 31.1 33.2 8.1 12:9 13,5 155.4 0 3.4 
Albercas 236 166 225 43 97/22 10.1 1838 195) 182.3 0 Ox 
Yerbab 324 228 22 43.3 38.6 10.7 13.6 15.6 183.1 0 9.3 
Urepetiro 386 270 29.4 47.5 SIE) 92 13.6 DON) 16078 0 0) 
Bosque 196 138 20.4 40.6 40.3 8.2 14.7 12.4 195.0 0 1.1 
Calicanto 232 163 25 20.3 24 9.6 15.6 12.8 1295) 0 Sl 
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Cuadro 2. Resultados de análisis físico-químicos a muestras de agua de pozos. 


e o a E K Mg Ca cr | HCO; | CO- | sos 
(ppm) (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) 
Morelos 120.7 84.7 17.5 20.1 15.2 6.2 8.4 Lal 101.7 0 2.2 
Bombas 404 283 197 38.5 36.4 10.4 2DAS 14.2 182.6 0 20 
Llano 246 173 20 31.9 28.8 9.9 14.3 19.1 138.2 0 8.9 
Cuates 144 100.9 19 32.6 812 92 122 1575) 134.0 0 8.6 
Palmar 251 176 21 32.9 30.4 9:5 12.6 19.9 137.1 0 9 
Tlazazalca 363 255 23.5 46 Sal 10.9 11512 15.6 211.0 0 pia 
Parocuello 428 300 27.5 48.2 52.6 9.2 15 37.6 189.6 0 7A 
Tenería 207 147 ZO 40.6 50.8 8.2 11545) PANES) 188.8 0 2.4 
Camelinas 365 254 22.6 539:2 56.2 10.1 11.3 53.9 182.9 0 5.9 
Loma 389 273 21.6 S28) 45.3 1187 16.8 11872 199.1 0 PS) 
Duarte 490 343 20.2 40.1 51,3 19.2 14.3 10.3 225.0 0 20.5 
Moral 492 345 2 35.8 5215 17.4 289) 17.0 210.0 0 26.8 
Maguey 530 370 19.9 34.5 45.5 18.3 27.6 17.7 235.0 0 5 
Tosca 405 286 24.5 66.4 SS) 10.1 19.8 45.4 205.0 0 17S 
Llano 512 359 20.7 45.3 35.6 18.3 25.9 18.4 245.0 0 19.2 
Cofradía 145 102 ALL 41.7 20.7 8.2 14.2 10.6 165.0 0 2.4 
Porvenir 414 293 29.6 55.3 70.2 4.8 16.6 18.5 235.0 0 6.8 
Floresta S9l 212 27 45.8 50.6 ES 16.5 EZ LS 0 LA 
DUM 557 266 23.5 45.5 50.2 16.5 42.3 21.3 310.0 0 10.6 
Gortari 766 452 2242) 58 95 2L2L E) DES 43.2 369.6 0 40.6 
Arbole 520 333 23.2 50.2 72 22.9 45.5 22.0 324.0 0 39.2 
Jericó 366 235 26.7 47.8 DOLO, 910 37.6 Ez 268.2 0 DOES) 
Villa 870 609 31.9 55.9 61.5 17.5 26.1 21.6 298.0 0 50 
Agustín AOS 774 2 O) 58.3 2L2LE) 38.5 DS) 294.0 0 SI 
Riveras 741 518 23.1 52.7 54.5 20.1 36.9 29.5 295.1 0 25.3 
Peñita 1 021 ALS 25.4 61.4 65.9 28.6 38.6 40.8 354.5 0 15.9 
Realejo 717 503 24.5 52.3 58.3 17.4 45.4 22.3 315.4 0 19 
Nuevo 476 334 30.3 45 55 12.8 29 9.6 280.3 0 0 
Unidad 141 99 22.7 28.7 23.8 6.4 8.5 5.7 130.4 0 0 
En Bosque 9/2 261 24.5 34.3 36.9 8.6 12O) 13.8 164.9 0 2.4 
3 Entierro 203 142 22.3 28.2 32.5 8.2 10 12.8 139.9 0 1.4 
3 Pivac 1 417 997 20.6 52 81.3 42 55.6 26.6 250.9 0 191 
E Hacienda 305 240 27.2 42.1 45.6 72 12.2 8.2 208.9 0 0 
E Calera 570 399 SS) 45 65 9.5 16.8 10.3 244.0 0 0 
E Luis G. 1520 1 367 22.9 330 84.5 23.6 63.4 58.8 887.1 0 125 
E Ixtlán 830 580 43.5 1549) 48.2 2d) 14105 74.4 363.7 0 32.8 
8 Chinela 2 030 1510 23.3 380.8 95.2 26.6 65.1 355.3 450.0 0 149 
E Conce 750 523 34.2 190.7 66.4 ES 13.4 65.4 425.0 0 105) 
= Geyser 2 630 1860 43 83.7 84 11 18.5 25.6 330.8 0 7.9 
E Sauceda 368 255 325 52.6 2 6.4 1221 ez 210.3 8.9 2.6 
A Sauceda 2 609 426 34.4 619.7 166.7 3.7 21.3 723.9 67.5 17.6 348 
3 Rastro 1 006 706 24.6 212.4 2087) TAS 38.4 26.9 639.1 0 4 
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Cuadro 2 (continuación). Resultados de análisis físico-químicos a muestras de agua de pozos. 


O E O AS E K Mg Ca Cr | HCO; | CO | SO4 
(ppm) (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) | (mg/L) 
Esquinas 436 308 26 78.2 7.8 4.9 10.0 31.9 201.3 0 9.6 
J. Rosario 673 473 25 115 15.6 97 36.1 32 482.0 0 14.4 
L. Ordaz 1315 935 23.5 133.4 11:7 45.0 140.3 53.2 11226 0 19.2 
Tule 1530 1145 DSi 207 117 45.0 128.3 39.0 1 134.8 0 14.4 
Luis Ríos 1 444 1016 22.3 186.3 11.7 36.5 120.2 39.0 933.5 0 24.0 
J. Zepeda Z. 1510 1189 211€) 243.8 1.7 36.5 150.3 42.6 1 104.3 0 187.22 
A. Gamez 1500 1079 23.3 225.4 11.7 36.5 108.2 81.6 982.3 0 9.6 
Pajacu 1412 991 24.4 151.8 101,7 42.6 156.3 032 OAES 0 14.4 
M. Estrada 592 416 24.3 92 7.8 9.7 36.1 28.4 384.4 0 14.4 
E. Barajas 1 457 1030 2 230 1.7 34.0 24.0 46.1 9PES 0 28.8 
A. Rodríguez 1377 973 24 167.9 11.7 40.1 112.2 53.2 890.7 0 9.6 
M. Navarro 1520 1 080 25 209.3 161.7 38.9 16622 532 1 018.9 0 9.6 
Gregorio 315 230 24 48.3 7.8 8.5 20.0 24.8 183.0 0 9.6 
J. Morales iS) 827 21.9 140.3 1.7 32.8 12.0 SIR9 616.2 0 38.4 
R. Garibay 885 625 23.5 119.6 7.8 19.5 88.2 31.9 347.8 0 172.9 
Pajacu 345 241 3212 43.7 7.8 97 36.1 14.2 29) 0 4.8 
Viejo 367 257 32.7 55.2 7.8 21.9 12.0 21.3 231.8 0 4.8 
Morales 875 616 26.8 121.9 7.8 2 12.0 SS) 488.1 0 33.6 
R. Ríos 475 330 36.4 71.3 11.7 4.9 16.0 70.9 207.4 0 9.6 
A. Ríos 495 345 34.5 ES 1.7 6.1 10.0 63.8 183.0 0 24.0 
Higuera 475 330 33.9 69 11.7 6.1 12.0 53.2 207.4 0 19.2 
G. Morales 455 Si 972 69 7.8 6.1 10.0 5512 207.4 0 9.6 
J. Padilla 512 359 32 71.3 11.7 8.5 16.0 49.6 237.9 0 28.8 
E. Jiménez 360 202 SS) 50.6 7.8 TIM 12.0 24.8 201.6 0 4.8 
Paracho 339 239 30.2 48.3 7.8 8.5 12.0 17.7 158.6 0 4.8 
M. Valdez 1 580 1 ES DOS 354.2 7.8 24.3 78.2 46.1 META 0 14.4 
A. Bravo 1810 1392 23.9 324.3 11.7 34.0 126.3 53.2 1 336.1 0 115.3 
R. Toribio 1850 1 418 24.7 301.3 16,7 34.0 124.2 166.7 1 049.4 8 (557 
R. Ávalos S. 1520 1200 30.6 239.2 11:7 30.4 116.2 102.8 988.4 0 57.6 
Meza 935 656 229 119.6 7.8 20.7 64.1 53.2 555.2 3 48.0 
J. González 1540 1115 25:3 216.2 11.7 32.8 114.2 60.3 1 098.2 0 28.8 
Buenrostro (192 835 24.7 147.2 7 26.8 92.2 63.8 780.9 8 14.4 
V. Arévalo 1540 1129 25.3 188.6 11.7 37.7 146.3 46.1 11775 3 9.6 
D. Gómez 1 500 1 054 DO) 200.1 7.8 38.9 102.2 5 1 000.6 0 67.2 


Se obtuvieron rangos de valores de CE, 
de 101 a 2 630 uS/cm; los valores más bajos 
que representan una mejor calidad química 
del agua subterránea se distribuyen en las 
estribaciones de las sierras de mayor altitud, 
localizados en la porción sur-sureste de la 
cuenca; en tanto que los valores más altos se 


registraron al NW del área de estudio, que se 
relacionan con la zona de influencia del campo 
geotérmico de Ixtlán de los Hervores, así como 
con la porción central orientada en el sentido 
de la falla geológica Pajacuarán. 

En relación con STD, los rangos oscilaron 
entre 60 y 1 800 mg+1*, coincidiendo los valores 
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Figura 3. Mapas de isovalores para los parámetros de CE, STD y T. 


más bajos con las estribaciones de los cerros de 
mayor altitud, aumentando su concentración 
de manera gradual hacia los valles; lo anterior 
es una guía de la presencia de agua de mejor 
calidad asociada con el tipo de acuífero en 
explotación. 

El área de mayor concentración se relaciona 
con la zona de influencia del campo geotérmico 
de Ixtlán de los hervores, así como con la 
porción central orientada en el sentido de la 
falla geológica Pajacuarán. 

La relación de conductividad eléctrica (CE) 
es directamente proporcional al contenido de 
sólidos totales disueltos en el agua subterránea 
y éstos se incrementan con la evolución del agua 
subterránea, o por efecto de la contaminación 
natural o antropogénica. 

La temperatura puede 
incremento o decremento de la solubilidad de 


influir en el 


las sustancias presentes (Porras-Martín et al., 


1985), en específico su conductividad eléctrica; 
en la cuenca, los valores obtenidos fueron 
entre 17.5 y 43 *C. Las temperaturas más frías 
se ubican hacia las estribaciones volcánicas, 
coincidiendo con la alta concentración de 
manantiales. Los valores mayores se asocian 
con el sistema estructural de la porción norte, 
que dan lugar a condiciones geotérmicas 
locales. En esencia, de 
temperatura están relacionadas con el esquema 


estas anomalías 
geológico-estructural de la región. 

La figura 4 muestra la distribución de 
isovalores de RAS para la cuenca. Los valores 
obtenidos de RAS oscilaron entre 1 y 32 meq'1*, 
lo que refleja que, en términos del elemento 
sodio, la mayoría de los aprovechamientos 
representan un grado bajo de peligrosidad. 
Destaca la zona de la ciénaga con valores 
máximos y la parte alta de la cuenca con los 
valores menores. 
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Figura 4. Mapa de isovalores de la relación de adsorción de sodio (RAS) en la cuenca. 


En relación con la clase de agua presente, la 
figura 5 muestra la ubicación de manantiales 
y pozos, respectivamente, en el diagrama 
referido. 

El agua 
principalmente como del tipo C1-SI, con bajo 


de manantiales se clasifica 
peligro de salinización y sodización, aptas para 
utilizarse en el riego agrícola sin restricción. 
Es decir, son aguas de excelente calidad 


química. Las clases de agua para los pozos 


son predominantemente del tipo (2-51 con 
salinidad media y bajas en sodio, pudiéndose 
utilizar en mayor medida en riego agrícola. 
La mayoría de ellos se ubican en un ambiente 
geológico de tipo fluvial. 

En la región lacustre se encontraron aguas 
del tipo C3-S1, con peligrosidad alta a la 
salinización, por lo cual su uso debe restringirse 
al riego de productos tolerantes a dicha 
condición. Finalmente, el ambiente lacustre 
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Figura 5. Ubicación de los sitios de muestreo de aguas subterráneas en el diagrama de H. Riverside: a) manantiales y b) pozos. 
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influenciado por hidrotermalismo, la clase de 
agua predominante, es (3-53, con salinidad 
muy alta y peligro de sodio alto, por lo que su 
uso en la agricultura se debe restringir. 

Lo anterior demuestra la interacción 
agua-roca como mecanismo de control de la 


calidad química del agua subterránea en la 


cuenca hidrológica del río Duero para uso 
primordialmente agrícola (figura 6). 

Por otra parte y con la utilización del 
diagrama de Piper, se puntualizaron las 
distintas facies químicas de agua presentes en 
la cuenca. La figura 7 muestra su distribución 
en el diagrama respectivo. 
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Figura 6. Distribución de las distintas clases de agua en la cuenca asociado con el ambiente geológico. 
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Figura 7. Ubicación del muestreo hidrogeoquímico en el diagrama de Piper: a) manantiales y b) pozos. 
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En general, tanto el agua de manantiales 
como de pozos se clasifica como de tipo 
bicarbonatada-sódica, que corresponde con 
aguas poco evolucionadas, consideradas como 
típicas de ambientes geológicos volcánicos. 


Agua superficial 


La figura 8 muestra la ubicación de los sitios 
de muestreo a lo largo del cauce del río Duero. 

En cuanto a las variables ambientales, 
los valores de oxígeno disuelto fueron 
relativamente bajos a lo largo de todo el río 
(cuadro 3). El pH resultó tendiente a la acidez, 
aunque con cantidades muy aproximadas al 
valor neutro. Los valores de la conductividad 
eléctrica aumentaron de manera gradual, 
conforme se dirige de la parte alta a la parte 
baja de la cuenca. 

De manera similar, la alcalinidad cambia de 
aguas blandas en la parte superior (0-75 mg;1*) 
a aguas duras en la parte baja de la cuenca 
(150-300 mg1*, con la clasificación de Sawyer 


y McCarty, 1967). Otras variables, como el 
DBO,, presentan patrones de fluctuaciones, los 
cuales muestran procesos de dilución hacia la 
parte alta y procesos de sedimentación en la 
parte baja. Por su parte, aquellas variables que 
representan nutrientes derivados de drenajes 
y de residuos de fertilizantes de la agricultura 
indican sitios de mayor impacto por estas 
actividades, en especial en los valles (cuadro 3). 

La interpretación de los valores ICA en los 
ocho manantiales analizados, en términos de 
su uso en la categoría de agua potable, indica 
que sólo el manantial de Carapan (origen del 
río Duero) muestra condiciones de excelente 
calidad (92 puntos; figura 9). La localidad del 
manantial de La Luz presenta contaminación 
leve. El resto de los manantiales, como es 
el caso de Camécuaro, Orandino y San José 
de Gracia, se clasifican como contaminados, 
destacando en particular la presencia de 
bacterias coliformes fecales. 

En lo que respecta al cauce del río, se 
implementó la interpretación del índice en el 
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Figura 8. Localización de los sitios de muestreo en el cauce del río Duero. 
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Cuadro 3. Variables ambientales del río Duero: Temp = temperatura, OD = oxígeno disuelto, STS = sólidos totales suspendidos, 
Alc = alcalinidad, DBO = demanda bioquímica de oxígeno a cinco días, DQO = demanda química de oxígeno, 
C. fecales = coliformes fecales e ICA = índice de calidad del agua. 


Alc total 
auesHa a Se Pa Ea pa a dal a pote ps on IEA 
CaCO, 
Manantiales 
Carapan 1618 |862| 51 | 01 |644| 1783 0 0 0 0 0 92 
Tacuro 16.6 |573| 408 | 0.11 |692| 1898 | 16 | 24 0.7 7.97 24 000 60 
San José de Gracia | 17.97 |75.5 | 5.13 | 0.08 |6.45 | 58.2 0 0 0.34 | 1683 930 68 
Chilchota 16.73 | 57 | 402 | o131 | 67 | 1791 | 11 | 16 | 079 | 1905 240 000 56 
Santo Tomás 163 (555 | 4 012 |6.65| 184,5 0 0 0.75 3.54 2.300 64 
Camécuaro 2033 |618 | 41 | 013 [693 | 852 | 11 7 0.37 93 2.400 64 
Orandino A E O O e o 0 0 0.32 3.99 230 69 
La Luz 2045 | 66.6 | 437 | 0.11 |6.78| 82 0 0 0.34 3.99 230 71 
Cauce del río 
Presa La Yerbabuena | 23.92 | 46.7 | 2.91 12.08 |7.17 | 146.7 32 35 1.14 0 40 63 
Río Tlazazalca 1992 |394 | 262 | 251 |689| 852 | 246 | 71 3.5 17.72 460 000 38 
Presa Urepetiro 2338 |33.7 | 212 | 1533 |6.95| 1456 | 69 | 45 0.69 7.53 4.300 47 
Duero ADV Gpe. | 1759 |54.9 | 3.84 | 116 |681 | 9 43 6 0.91 7.09 15 000 55 
Duero DV Gpe. 1879 | 63 | 432 | 1183 [659 | 94 43 | 62 | 081 7.97 23 000 55 
Río Sagrario 19.87 | 39.4 | 254 | 18.07 |6.66 | 1165 | 94 | 37 3.64 | 28.35 2400000 | 39 
Río Urepetiro 192 |623| 424 | 1899 | 66 | 1009 | 86 | 36 0.82 4.43 9 000 55 
Las Adjuntas 1788 |57.2| 403 | 1241 (663 | 99.2 | 39 13 0.91 6.09 93 000 54 
Unión Tangacícuaro | 17.92 | 58 | 401 | 1228 |6.48 | 938 | 36 | 98 | 103 6.65 23 000 54 
El Platanal 19.16 | 67.1| 457 | 14.74 (682 | 1738 | 35 15 0.98 | 1595 43 000 56 
e 1857 1685 | 472 | 1222 [6.71 | 1025 | 05 | 95 | 0.75 14.62 9 000 62 
pera Dueto ome esy | 427 ss les | eso | a 32 0.8 15.51 23 000 53 
Río Duero 19.08 | 52 | 3.55 | 12.62 |6.57 | 95 45 | 205 | 0.91 8.64 43 000 51 
Río Celio 2026 |605 | 3.98 | 1131 |6.78| 85 08 | 13 | 059 5.98 43 000 60 
Puente Los Espinos | 20.78 | 76.1 | 5.02 | 13.54 | 6.58 | 171.3 10.8 26.5 4.78 7.53 93 000 51 
a eii 2053 | 69.2 | 4.55 | 13.58 [6.51 | 1569 | 44 | 23 17 15.95 93 000 52 
3 
3 pan cen 25.7 |84.8| 5.04 | 39.04 |7.56 | 280.5 | 535 | 122 | 496 | 14.18 24000 48 
E Dren A 2081 | 65.7 | 416 | 15.02 | 6.47 | 1254 6 10.8 | 196 7.97 93 000 51 
E Dren Colongo 2085 |51.1| 3.37 | 14.58 |6.61 | 1079 | 57 | 25 1.16 7.53 460 000 49 
E La Estanzuela 2093 |463 | 299 | 1473 | 66 | 1441 | 69 | 20 | 202 | 1329 93 000 45 
3 Dren Camucuato | 2321 | 14 | 08 | 4236 |6.78| 1144 4 40 2.28 3.54 110 000 44 
E San Cristóbal 2157 [639 | 206 | 161 [68 | 1572 | 52 | 26 2.59 8.2 110 000 51 
E El Capulín 22.26 |38.4 | 246 | 1644 |6.65 | 1492 | 84 | 55 1.39 731 46 000 44 
3 Puente Cumuato 2453 |55.7 | 332 | 2397 |6.81 | 1305 | 62 | 35.5 | 2.0175 | 653 230 53 
- Duero-Ibarra 23.66 |59.9 | 3.68 | 182 |659| 136 32 | 175 | 0.89 3.32 40 64 
S Río Lerma 2342 |a7a | 248 | 22.17 [6.66 | 2015 | 39 | 268 | 3.04 8.64 110 000 46 
=> Unión Duero-Lerma | 23.43 | 25.4 | 153 | 21.47 |6.72 | 106 59 | 215 | 139 | 10.19 46 000 42 
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Figura 9. Variación del ICA para el uso de agua potable en los manantiales del río Duero. 


aspecto agrícola, quedando los diferentes sitios 
en la condición de levemente contaminado 
y contaminado. Destacan como los sitios de 
mayor contaminación el río Tlazazalca y río 
El Sagrario en Tangancícuaro, este último reci- 
biendo las descargas domésticas de la capital 
municipal y de actividades agropecuarias, lo 
que degrada su calidad del agua de manera 
importante. Otros lugares con valores bajos del 
ICA son la presa Urepetiro, La Estanzuela, dren 
Camucuato y El Capulín (figura 10). 

El análisis de las fluctuaciones del valor 
del ICA a lo largo del trayecto del río refleja 
procesos de concentración de contaminantes y 
de su dilución por la entrada de aguas menos 
deterioradas, como de los manantiales. Cabe 
destacar que en la parte final del trayecto, antes 
de la desembocadura en el río Lerma, cambia 
de contaminado en El Capulín, a levemente 
contaminado en Cumuato e Ibarra, porque esta 
zona es más léntica y se tiene un proceso de 
sedimentación. 


Un comparativo realizado para los últimos 
20 años muestra el avance del deterioro que 
manifiesta el río Duero en términos de calidad 
de agua. Por ejemplo, el oxígeno disuelto en 
algunos tramos distribuidos a lo largo del río 
disminuye, rebasando los límites de tolerancia 
de los organismos acuáticos (La Estanzuela y 
Briseñas). 
tal 
comparándose valores actuales con datos 
reportados por (1988) para esas 
localidades. También se registra un incremento 


La figura 11 muestra tendencia, 


López 


en los valores de turbidez de hasta dos órdenes 
de magnitud en algunos sitios (López obs. 
pers.), lo que está directamente relacionado con 
la erosión y descargas. 

Por su parte, los nutrientes también 
aumentaron al comparar cuatro años. Destaca 
en particular que la diferencia entre los años 
se da principalmente entre los tres primeros, 
que fueron monitoreados de manera seguida 


en la década de 1980 y el valor obtenido en 
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Figura 10. Variación del ICA para el uso de agua de riego en diferentes sitios del cauce del río Duero. 
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el 2009, indicando una diferencia estadística 
importante (F = 15.73, p = 0.0001) (cuadro 4). 


Discusión 


Los problemas de contaminación de cuencas en 
México se han agudizado al paso del tiempo. 
Su deterioro se ve reflejado en los afluentes 
y reservorios, los cuales indican el estado 
de salud de la cuenca, de lo que pasa en ella. 
Esto es especialmente importante debido a 
que gran parte de la actividad económica 
que se desarrolla depende de la calidad del 
agua disponible en las cuencas. Para describir 
esta calidad se pueden tener diferentes 
aproximaciones, desde la clasificación en 
escalas de uno hasta tres parámetros, el análisis 
de diferentes parámetros de una manera 
individual o multivariada, por medio de 
estadísticos de agrupación u ordenación, hasta 
la implementación de índices. 

El organismo encargado de monitorear la 
calidad del agua en México es la Conagua, la 
cual emplea índices de acuerdo con parámetros 
químicos. Antes del año 2002, la Conagua 
utilizó un índice de calidad del agua basado 
en un promedio ponderado de los índices de 
calidad de 18 parámetros representados en una 
escala de O a 100. El uso de este método colocó 


73% de las corrientes de agua superficial en la 
condición de contaminado (Conagua, 2001). 
Posterior al 2002, la Conagua decidió cam- 
biar la forma en que se establece la calidad del 
agua de aguas superficiales. Para ello utiliza 
principalmente de manera independiente los 
parámetros de demanda bioquímica de oxígeno 
(DBO), demanda química de oxígeno (DQO) y 
sólidos suspendidos totales (SST) (Conagua, 
2011). El valor resultante de estos parámetros 
se contrasta con un criterio de clasificación de la 
calidad del agua. A pesar de que es una forma 
más sencilla, rápida y efectiva —en términos 
de costos— de obtener una valoración de dicha 
calidad, también presenta ciertas limitantes. 
La principal limitante en el caso del DBO y 
DQO es que determina aquellos efectos debido 
a los cambios en las descargas orgánicas, 
implicando que la calidad está siendo definida 
en relación con la contaminación por materia 
orgánica. Los SST indican, sin discriminar, 
la presencia de diferentes elementos: limo, 
arcillas, sedimentos de escorrentías, plancton, 
desechos industriales y de aguas residuales, 
entre otros, presentes en el agua. Esto crea una 
tendencia a mostrar que el panorama en general 
es mejor de lo que realmente está (SEPA, 2013). 
Para corroborar lo anterior, cuando se realiza el 
análisis de la situación de la calidad del agua 


Cuadro 4. Análisis Tuckey HSD para comparaciones múltiples en cuatro años de la variable nitratos en seis sitios de muestreo. 


(D) Años () Años (-D Error estándar Sig. Límite inferior | Límite superior 
1985 0.2298 0.4185 0.946 -0.9416 1.4013 
1984 1986 -0.4408 0.4185 0.721 -1.6323 0.7306 
2009 -1.7632* 0.4185 0.002 -2.9346 -0.5917 
1984 -0.2298 0.4185 0.946 -1.4013 0.9416 
1985 1986 -6.71E-01 0.4185 0.4 -1.8421 0.5008 
2009 -1.9930* 0.4185 0.001 -3.1644 -0.8216 
1984 0.4408 0.4185 0.721 -0.7306 1.6123 
1986 1985 6.71E-01 0.4185 0.4 -0.5008 1.8421 
2009 -1.3223* 0.4185 0.23 -2.4938 -0.1509 
1984 1.7632* 0.4185 0.002 0.5917 2.9346 
ld 1985 1.9930* 0.4185 0.001 0.8216 3.1644 
1986 LEPE 0.4185 0.023 0.1509 2.4938 


* significativo al 0.05 
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con estos criterios en el año 2004, se tuvo un 
cambio a tan sólo un 17.4% de los ecosistemas 
acuáticos con algún grado de contaminación 
(Conagua, 2004). 

En el caso del río Duero, también 
se identifican diferencias al aplicar los 
distintos criterios con respecto al ICA. En 
los manantiales, el DBO, DQO y SST dan la 
clasificación de Excelente, mientras que en el 
ICA sólo el de Carapan alcanza dicha categoría, 
y la mayoría está al nivel de Contaminado. Esto 
se relaciona en buena medida con la presencia 
de coliformes fecales, incluso Escherichia coli, 
que se utilizan como indicadores de posible 
contaminación por patógenos. Además, los 
nutrientes —nitratos y fosfatos— en algunos 
casos alcanzaron valores elevados y están 
ligados con procesos de eutroficación. 

En el cauce del río, aun con un criterio 
menos estricto —uso agrícola— y equiparando 
las categorías  —dándoles el nombre 
equivalente—, no hay concordancias entre 
los indicadores. Incluso se llega en el sitio de 
la unión del Duero con el Lerma a tener cada 
indicador una categoría diferentes (DBO: Buena 
Calidad; DQO: Aceptable; SST: Excelente, e 
ICA: Contaminada). También cambian las 
proporciones por criterio; mientras que en el 
DBO se tiene un 26% de sitios catalogado como 
aceptable, en el DQO es 48% y en el ICA 63%. 
Por su parte, 4% de sitios está contaminado de 
acuerdo con el DBO; 19% en el DQO, y 37% en 
el ICA. 

Esto no significa que se esté siguiendo 
un Criterio erróneo, de antemano se ha 
reconocido el valor del DBO como indicador 
y se ha registrado durante extensos periodos 
de tiempo en diferentes países como el caso 
de Inglaterra (Callow, 1990). Ellos incluyen 
además parámetros físicos con los que se 
cuentan en México, como el flujo del río y la 
temperatura, al igual que nutrientes como el 
nitrógeno amoniacal y el fósforo total. Ambos 
se relacionan con fertilizantes, ganadería, y 
aguas residuales industriales y domésticas 
(SEPA, 2013). Por su parte, si bien los ICA 
proveen valores significativos que indican la 
salud del recurso hídrico, y crean un criterio 


para medir y calcular la calidad del agua, 
también presentan imprecisiones. 

El ICA propuesto (NSFWQI) se ha señalado 
como no muy objetivo, porque un grupo de 
expertos califica y escala los parámetros de 
calidad del agua que se utilizan —apodado 
como método Delphi— y dicho escalamiento 
puede cambiar en la apreciación de otro grupo 
(Lumb, 2011). Incluso el índice se ha modifica- 
do para atender las necesidades específicas 
de los estados, como en el caso de Oregon 
(OWQI; Cude, 2008) e lowa en los Estados 
Unidos (IDNR, 2013). Su implementación en el 
presente trabajo parte no sólo de su aún común 
implementación, además de su formulación 
multiplicativa que le da robustez en cuanto al 
mejor ajuste a las opiniones de los expertos, 
sino a la disponibilidad de la información para 
incorporar un mayor número de sitios. 

Más aún, el ICA en sus diferentes moda- 
lidades proporciona una visión parcial sobre 
la salud de los ecosistemas acuáticos. Cuando 
se desea evaluar el estado de los ecosistemas 
es importante incorporar de manera adicional 
el uso de índices biológicos; ente los más 
utilizados se encuentran el Índice Biótico 
Estendido (IBE) y el Índice de Integridad 
Biótica (IBID), basado en macroinvertebrados 
y peces, respectivamente. Estos índices se 
han propuesto en diferentes ocasiones para 
cuerpos de agua mexicanos (Weigel et al., 2002; 
López et al., 2007; Mathuriau et al., 2011). Ante 
esta situación se debe establecer un esfuerzo 
conjunto para monitorear de forma integral los 
ecosistemas acuáticos, aunque no sea de manera 
continua, pero sí con cierta periodicidad. 

En la cuenca del río Duero es muy clara 
la tendencia a un creciente proceso de de- 
gradación, según se demostró en el análisis 
temporal. Esto es crítico, dado que la principal 
actividad económica de la región es la 
agricultura. En el ámbito nacional, se distingue 
por ser la principal productora de fresa, y en 
los últimos tres años, en cultivos de frutillas, 
como zarzamora, arándanos y frambuesa. 

Entre las problemáticas que afectan a estos 
cultivos se tiene principalmente la calidad 
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del agua para su riego. El mercado externo 
requiere productos con alta inocuidad. Tiempo 
atrás, el agua de riego utilizada era la del río, 
sin embargo hoy en día esto es imposible por el 
deterioro de su calidad. 

La problemática sentida en la calidad del 
agua superficial fue el detonador de la apertura 
de una gran cantidad de pozos profundos para 
explotar el agua subterránea, pasando de un 
registro de 150 aprovechamientos para el año 
2000 a 805 en la actualidad (REPDA-Conagua). 
La apertura de estas concesiones se ha realiza- 
do sin planeación alguna, como lo demuestra 
la alta concentración de pozos, pudiéndose 
encontrar hasta 20 pozos por áreas de 4 km?. 
La figura 12 muestra esta problemática, 
definiéndose tres niveles: baja (0-2), moderada 
(3-4) y alta (> 4); como consecuencia, se inicia 
un círculo vicioso de perforación de pozos que 
puede propiciar su agotamiento. 


Aunado a lo anterior está la calidad del agua 
subterránea, que si bien en términos generales 
es buena para la cuenca, hay sitios donde la 
gran presencia de sales (en los alrededores de 
Ixtlán) llega a causar problemas de salinización 
de los suelos. En estos casos se requiere un pre- 
tratamiento. 


Conclusiones 


El estudio hidrogeoquímico del agua subte- 
rránea en la cuenca demostró la interacción 
agua-roca como mecanismo de control de su 
calidad química para uso primordialmente 
agrícola. La clase de agua predominante es C1- 
S1, baja en salinidad y en sodio, pudiéndose 
utilizar en la mayor parte de los cultivos, así 
como en cualquier tipo de suelos, sin desarrollar 
peligro por salinidad y sodicidad. 

La problemática principal en la cuenca con 
respecto al recurso hídrico es la mala calidad 


102% 20' 00” W 


102% 30' 00” W 


N 200 ¿OL ¿07 


N 00 ,0 ¿07 


KO] Cuenca del río Duero 
*X Cabeceras municipales 
Pozos 
UC Areas urbanas 
—— Río Duero 
MM Cuerpos de agua 
O Límite estatal 
Límites municipales 
Densidad de pozos 
Cuadrante de 2x 2 km 
[Adecuado (1-2 pozos) 
| [A Medianamente adecuado (3-4 pozos, 
IE No adecuado (5-24 pozos) 


N 200 ¿05 ¿61 


102% 20' 00” W 


102% 10' 00” W 1022 00' 00” W 1019 50' 00” W 


N ¿,00 01 ¿07 


N 00.0 ¿07 


N 200 05 ¿61 


102% 00' 00” W 


102% 10' 00” W 101% 50' 00” W 


Figura 12. Densidad de pozos por cuadrantes de 2 x 2 km en la cuenca del río Duero. La imagen muestra la marcada sobre 


posición de dominios entre ellos. 
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del agua del río Duero, relacionada con las 
descargas directas de aguas residuales, por lo 
general de origen urbano, detectadas a lo largo 
del río. Es resultado de la falta de compro- 
miso por parte de las autoridades munici- 
pales por cumplir con su responsabilidad de 
saneamiento y se complica de acuerdo con el 
esquema normativo, que permite se mezcle 
el agua de manantiales con las descargas 
urbanas. También un factor decisivo es la 
poca coordinación entre autoridades en los 
tres niveles de gobierno para el desarrollo de 
acciones, y una escasa o nula cultura del agua 
entre los pobladores de la cuenca. 
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Abstract 


SILVA, J.T, MONCAYO, R., OCHOA, S., ESTRADA, F, CRUZ-CÁRDENAS, G., 
ESCALERA, C., VILLALPANDO, E, NAVA, J]., RAMOS, A. €; LÓPEZ, M. Groundwater 
and surface water quality in the Duero River basin, Michoacan. Water Technology and 
Sciences (in Spanish). Vol. IV, No. 5, November-December, 2013, pp. 127-146. 


This work describes the chemical quality of groundwater and surface water resources in the 
Duero River basin. The groundwater was analyzed by a hydrogeochemical study of 97 water 
sources (20 springs and 76 wells) during the dry season. The Water Quality Index (NSFWQD 
was used for 35 surface water sites, with eight variables divided according to springs and the 
river. The overall groundwater quality was good and was associated with the basin's rocks and 
geology. The predominant class was C1-S1, with low salinity and sodium, useful for most crops 
and any soil types, with no risks due to salinity and sodicity. Based on the NSFWQI values 
obtained for potable water, most of the springs were classified as contaminated, except for the 
Carapan spring (excellent quality, 92 points). On the main river channel, with agriculture 
use, 75% of the sites were determined to be contaminated and the rest slightly contaminated. 
The comparison of nitrates over 20 years presented a statistically significant difference (F = 
15.73, p =0.0001), indicating a significant increase in nutrients in the river. Finally, despite a 
positive water balance, the increasing deterioration in the quality of the water in the river leads 
to groundwater extraction and the resulting overexploitation of the aquifer. 


Keywords: hydrogeochemical analysis, pollution, overexploitation, water quality index. 
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